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Abstract 

We study the equilibrium statistical mechanics of classical two-dimensional 
Coulomb systems living on a pseudosphere (an infinite surface of constant 
negative curvature). The Coulomb potential created by one point charge 
exists and goes to zero at infinity. The pressure can be expanded as a series in 
integer powers of the density (the virial expansion). The correlation functions 
have a thermodynamic limit, and remarkably that limit is the same one for the 
Coulomb interaction and some other interaction law. However, special care is 
needed for defining a thermodynamic limit of the free energy density. There 
are sum rules expressing the property of perfect screening. These generic 
properties can be checked on the Debye-Hiickel approximation, and on two 
exactly solvable models : the one-component plasma and the two-component 
plasma, at some special temperature. 
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1 INTRODUCTION 



How the properties of a system are affected by the curvature of the space in 
which the system lives is a question which arises in general relativity. This is an 
incentive for studying simple models. 

The classical equilibrium statistical mechanics of Coulomb systems (charged par- 
ticles interacting by Coulomb's law) has been much studied. Coulomb systems have 
interesting generic properties, related to the screening effect. For two-dimensional 
systems (such as particles in a plane with a logarithmic interaction), there are exactly 
solvable models which are useful for checking and illustrating the generic properties. 
It might be of some interest to investigate what are the curvature effects. 

This problem has already been considered for Coulomb systems living on a sphere 
or an hypersphere 0, ||. One of the motivations for studying such systems 
is that they are convenient for numerical simulations, since they are finite without 
boundaries. However, that very finiteness has drawbacks from a fundamental point 
of view, since it is not possible to take the thermodynamic limit (infinite system) 
without suppressing the curvature at the same time ; thus one cannot get information 
about the effect of curvature on phenomena such as screening or phase transitions 
pertaining to infinite systems. 

In two dimensions, a surface of constant negative curvature may have the in- 
teresting feature of being both curved and infinite. The present paper is about 
Coulomb systems living on such a surface. In Section 2, salient properties of that 
surface are reviewed. The relevant Coulomb potential is discussed in Section 3. 
Section 4 is about general features of statistical mechanics on the curved surface. 
Screening effects are studied in Section 5. Then, more specific models are considered : 
the Debye-Huckel approximation in Section 6, the exactly solvable one-component 
plasma in Section 7, the exactly solvable two-component plasma in Section 8. 

2 PSEUDOSPHERE 

Let us recall a few properties of the surface of constant negative curvature called 
a pseudosphere |J. Such a surface is a two-dimensional manifold, the entirety of 
which cannot be embedded in three-dimensional Euclidean space. The surface can 
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be represented by a plane disk of radius 2a centered at the origin, the Poincare disk, 
with a metric such that 



ds> = dr2 + r2 y (2.1) 

in terms of the polar coordinates (r, y>). The geodesies are circle arcs normal to the 
boundary circle of the disk. The geodesic distance from any point inside the disk to 
the boundary circle is infinite : the size of the surface represented by the Poincare 
disk is infinite. The Gaussian curvature has the constant value — I /a 2 (our definition 
of the Gaussian curvature is such that it is 1/R 2 for a sphere of radius R). 
Another convenient set of coordinates is (r, tp) , with r defined through 

r = 2atanh^ (2.2) 
The corresponding expression for the metric is 

ds 2 = a 2 (dr 2 + sinh 2 r dp 2 ) (2.3) 
The geodesic distance between two points (r, <p) and (r', ip') is given by 

cosh — = cosh t cosh r' — sinh r sinh t cos(</2 — tp') (2-4) 

In particular, the geodesic distance to the origin is ar. 

This representation of a pseudosphere by a Poincare disk is very similar to the 
representation of a sphere of radius R by its stereographic projection from its North 
pole onto the plane tangent to its South pole 0. Here, R 2 is replaced by —a 2 , the 
spherical coordinate ir — 9 is replaced by r, trigonometric functions are replaced by 
hyperbolic ones. An important difference is that the (finite) sphere is represented 
by the whole plane, while the (infinite) surface of constant negative curvature is 
represented by the inside of the Poincare disk only. 

It should be noted that, for a system of particles living on a surface of constant 
negative curvature, having a uniform density n means that the average number of 
particles in the area element 
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is ndS and all points are equivalent. However, in the Poincare disk representation, 
there will be an apparent non-uniform density n[l — (r 2 /4a 2 )]" 2 which increases up 
to infinity as r approaches its upper value 2a. 

Through an appropriate conformal transformation, the Poincare disk can be 
mapped onto the upper half-plane, the Poincare half-plane. This provides another, 
often used, representation. Here, it will be more convenient to use the Poincare disk. 

3 THE COULOMB POTENTIAL 

The Laplacian is 



1 - 



4a 2 J 
1 d 
sinh t dr 



d 2 Id Id 2 
+ + 



Q r 2 r Q r r 2 Qp2 

B Id 2 
smh t— — h 



(3.1) 



dr ' sinh 2 r d<p 2 

The Coulomb potential v(s) at a geodesic distance s from a unit point charge is 
defined by 



Av(s) 



(3.2) 



together with the boundary condition that v(s) vanishes at infinity. In the present 
paper, unless otherwise specified, 5^ is the Dirac distribution on the curved mani- 
fold, i.e. such that 



6®(s)dS 



(3.3) 



where dS is an area element defined with the metric Q2.1|) . The solution of ( |3.2|) 
exists and is 



v(s) — — logtanh 



2a 



(3.4) 
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This result should be contrasted with the usual two-dimensional Coulomb poten- 
tial in a flat space, — \og(s/£), which cannot be made to vanish at infinity, whatever 
the choice of the constant length scale I might be. 

One should also note that, on a sphere, ( |3.2| ) has no solution, and one resorts 
to the choice |lj — log sin(s/2i?) which is the Coulomb potential created by a +1 
point charge and a —1 charge uniformly spread on the sphere, or alternatively the 
choice @ — logtan(s/2i?) which is the Coulomb potential created by a +1 point 
charge and a —1 point charge at the antipodal point. The difference is, of course, 
that, on a sphere, the field lines originating from a single point charge converge into 
another singularity at the antipodal point, while, on a pseudosphere, the field lines 
originating from a point charge spread away to infinity. 

4 STATISTICAL MECHANICS 

We now address the problem of building the classical equilibrium statistical me- 
chanics of an infinite system of charged particles living on a pseudosphere. Two par- 
ticles of charges and qj interact through the Coulomb potential —q%qj logtanh(sjj/2a), 
where Sy is the geodesic distance between these particles. We shall also consider the 
one-component plasma model, made of one species of positive particles of charge q 
and a uniformly negatively charged background. What are the generic properties to 
be expected ? 

4.1 Thermodynamic functions 

On a pseudosphere, there are problems with the definition of the thermodynamic 
functions. For instance, for defining the free energy density, the standard procedure 
would be to compute the free energy F in a domain of finite area S and to take the 
limit of F/S as the domain becomes infinite, keeping the particle number density 
n fixed. In the case of a flat system, the limit does not depend on the boundary 
conditions, because the boundary effects are of order S 1 ^ 2 while the bulk effects 
are of order S. However, in the present case of a pseudosphere, the length L of 
the boundary and the area S are of the same order. Indeed, a "geodesic disk" 
of radius s, i.e. the set of points which have a geodesic distance to some center 
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smaller than s, has an area S = Ana 2 sinh 2 (s/2a) while the length of its boundary is 
L = 27rasinh(s/a). As s — > oo, both S and L are of order exp(s/a). Therefore, the 
boundary contributions to F are important, and there is no unique thermodynamic 
limit of F/S. 

There are nevertheless reasonable ways of defining a free energy density or other 
thermodynamic functions. Some of these approaches will be described in the fol- 
lowing. Yet, it should be kept in mind that it is essential to get rid of the boundary 
effects. 



4.2 The virial expansion 

It is well-known that, for flat conducting Coulomb systems, the pressure cannot 
be expanded in integer powers of the density ; the pressure is not an analytical 
function of the density in a neighborhood of the origin. For those systems which, at 
low temperature, go into a dielectric phase through a Kosterlitz-Thouless transition, 
the non-analyticity can be used as a signature of the conducting phase || [|. 

It will now be shown that, on a pseudosphere, things appear to be different. 
There are strong indications that the virial expansion of the pressure in powers of 
the density exists for Coulomb systems which however will be shown in Sections 6 
and 7 to be conductors. 

At temperature T (as usual, we set (3 = 1/ksT with fee the Boltzmann constant), 
the virial expansion of the pressure p with respect to the total number density n 
would be 

oo 

(3p = n + Y J B k n k (4.1) 

k=2 

For simplicity, we shall only discuss in detail the second virial coefficient B2, for 
two models. 

4.2.1 Two-component plasma. This is a model made of two species of par- 
ticles of charges ±q, interacting through the Coulomb potential ( |3.4|) . The second 
virial coefficient is 



1 J ^-PiMs) + e P q Ms) _ 2 ] dS (4.2 
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where we have regrouped the Mayer bonds for like and unlike particles. The area 
element is 

dS = 2ixa sinh ^— \ ds (4-3) 

At large distances s, the square bracket in (|4~2| ) behaves like 4(/?g 2 ) 2 exp(— 2s/ a), 
and therefore the integral converges in spite of the sinh(s/a) in dS (some short-range 
potential should be added, if necessary, for making the integral convergent at short 
distances). This is to be contrasted with what happens for a flat system where B 2 
diverges. 

We conjecture that a similar analysis would show the higher-order virial coeffi- 
cients to be also finite. 

When (3q 2 < 2, no short-range potential is needed for making ( |4.2j ) conver- 
gent. With v(s) the pure Coulomb potential (|3.4|), B 2 can be explicitly computed. 
Through the change of variable t = tanh 2 (s/2a) and an integration by parts, one 
obtains 



Bo 



(3q 2 a 2 / : dt 

2 Jo 



i> 1 + 



l-t 



(4.4) 



2 J r \ 2 

where ip is the logarithmic derivative of the gamma function, which has an integral 
representation [10 which provides the second equality in (O). 



4.2.2 One-component plasma. This is a model made of one species of par- 
ticles of charge q and a uniform background charged with the opposite sign. The 
second virial coefficient now is 



Bo 



dS 



(4.5) 



--/ [e~^ - 1 + (3q 2 v(s) 

where we have added to the Mayer bond the contribution f3q 2 v(s) from the particle- 
background and background-background interactions. At large distances, the square 
bracket in ( f4.5| ) behaves like 2{(5q 2 ) 2 exp(— 2s/ a) and B 2 is finite. 

The explicit calculation of B 2 can be performed through the same change of 
variable as above t = tanh 2 (s/2a), with the result 
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Bo 



-1 1 _ ^9 2 /2 

-7l(3q 2 a 2 I — : : — dt 



where 7 = 0.577 • 



= —iij3q 2 a 2 
• is Euler's constant. 



l-t 



7 + V 1 + 



(3q 2 



(4.6) 



4.2.3 A possible definition of the thermodynamic functions. If indeed 
the virial expansion (|4.1|) exists, it provides, in principle, a definition of the pressure 
free from boundary effects. This happens because the thermodynamic limit of the 
virial coefficients Bk has been taken before the summation QO| ). The other ther- 
modynamic functions such as the free energy per particle / or the internal energy 
per particle u can be obtained through the relations 



p = n 



u 



(4.7) 



dn d/3 
Of course, in practice, the exact resummation of the expansion (|4. 1]) will not be 
feasible in general. 

The low-density expansion of the excess free energy per particle f exc starts as 



B 2 n + 



The corresponding expansions of the excess internal energy per particle are 
- for the two-component plasma 



n 2 2 
q a 



+ 



# 1 



(3q 2 



2 V 2 / - 

for the one-component plasma 



n + 



(4.9) 



Uexc = -7rq 2 a 2 



n + 



(4.10) 
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4.3 Correlations 



In contrast with the thermodynamic functions, the rz-body densities are ex- 
pected to have well-defined thermodynamic limits. This is related to the screening 
properties, to be discussed in Section 5 : screening makes the n-body densities inside 
the system, far away from the boundary, insensitive to the boundary conditions. 

Screening has a remarkable consequence : in the thermodynamic limit, the cor- 
relation functions are the same ones in a system with the Coulomb interaction 
— logtanh(s/2a) or when this interaction is replaced by — log sinh (s/2a). This can 
be shown as follows. 

iFiom ( |2.2j ) and (|2.4| ), it is found that, in terms of the complex coordinate 
z = rexp(iip) in the Poincare disk, the geodesic distance Sy between two particles 



located at Z{ and Zj is such that 



tanh — 
2a 



Zj )/2a 



1 - (ziZj/Aa 2 ) 



(4.11) 



and 



Si 



sinh ^ = N m (4.12) 

2a [1 - (r ?; /2a) 2 ] 1/2 [1 - (r,/2a) 2 ] 1/2 

For a system of particles of charges with the interaction — q^j logtanh(sy/2a), 

the total potential energy can be written 



-T^E^Og 

*¥=3 



Zj)/2a 



-gE^-log [l-(n/2a) 



ZiZj/Aa 2 

2 



{n/2af 



(4.13) 



The first two terms in the right-hand side of ( [4.13|) are the potential energy of 
a system of particles in a flat disk of radius 2a with ideal conductor walls at zero 
potential (the second-term is the interaction of each particle with its own image). On 
the other hand, for a system of particles with the interaction —qiqj logsinh(sjj/2a), 
the total potential energy can be written as 



Ho — — QiQn log — — — 

2^ MJ 2a 



$> 2 log \l-(r t /2af 



(4.14) 
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(the system has been assumed to be neutral, thus J2li — has been used). The 



first term in the right-hand side of ( 4.14 ) is the potential energy of a system of 
particles in a flat domain with plain hard walls. In both (OH) and (4.14), the last 
term is a same one-body potential associated with the curvature. Comparing ( |4.13| ) 



and ( 4.14 ), we see that the only difference between them is the nature of the walls. 
Since the apparent density in the Poincare disk (number of particles in d 2 z divided 
by d 2 z) becomes infinite at its boundary (because of the metric (|2.1| )), the apparent 
screening length will vanish at the boundary, and the correlation functions inside 
the Poincare disk will be independent of the nature of the wall. 
Similar consideration hold for a one-component plasma. 

In the following, for computing correlation functions, choosing the interaction 

— logtanh(s/2a) or — log sinh(s/2a) will just be a matter of convenience. In some 
cases, it will be explicitly checked that both interactions give the same correlation 
functions. The interaction — logsinh(s/2a) has the advantage of being the analog of 

— logsin(s/2i?) on a sphere |]], and that makes possible to use previously developed 
methods. 

Of course, both interactions have the same behavior — log(s/2a) in the flat sys- 
tem limit a — ► oo. 



5 SCREENING AND SUM RULES 

Screening is a characteristic property of conductors. Internal screening means 
that, at equilibrium, a particle of the system is surrounded by a polarization cloud 
of opposite charge. External screening means that, at equilibrium, an external 
charge introduced in the system is surrounded by a polarization cloud of opposite 
charge. For a flat system, internal screening and external screening of an infinites- 
imal charge, respectively, result in the two Stillinger-Lovett sum rules JTTJ] obeyed 
by the correlation functions. 

In this Section, assuming that screening holds, we derive the generalization of 
these Stillinger-Lovett sum rules to the case of a Coulomb system on a pseudosphere. 
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5.1 Internal screening 

Let us consider a system of particles of several species a. Each particle of species 
a carries a charge q a . The number density of species a is n a . The pair correlation 
functions are h a p(s) where s is the geodesic distance between the two particles. 
Internal screening is expressed by the sum rule 



fj2 n P IP h »/3(s) dS = -q a (5.1) 



5.2 External screening 

If an external point charge Q is introduced into the system, it induces a charge 
density Pq{s) at a distance s from this external charge. External screening means 
that 

J PQ (s)dS=-Q (5.2) 

In the case when Q is infinitesimal, linear response theory allows to transform 
(|5.2|) into a sum rule involving the correlation functions of the system in the absence 
of the external charge. Indeed, if the charge Q is at z, its interaction energy with 
the system is H int = Q<fi(z) where (f)(z) is the microscopic electric potential created 
at z by the system and the average induced charge density at z' is 

p Q (z') = -)3< p(z')H int >= -PQ < p(z')$(z) > (5.3) 
where p(z') is the microscopic charge density at z' . Using ( |5.2j ), one obtains the 



Carnie and Chan sum rule \T2 



(3 / < p(z')$(z) >dS' = 1 (5.4) 



Since < p{z')<j){z) > depends only on the distance s between z and z', ( |5.4j ) can be 
rewritten as 



2na 2 (3 I < p(0U(r) > sinhr dr = 1 (5.5) 
Jo 

and transformed by two successive integrations by parts into 



11 



Ana (5 / log cosh 
Jo V 2 



T 



A7Tr 3mhT Tr < ^ {T)> 



sinh rdr 



(5.6) 



where the Laplacian appears. From the Poisson equation, 



\ -rr— T- sinhr T- < ^(°)^( r ) >= ~ 27r < P(°)P( T ) > 
cr smhr ar ar 

Thus, we obtain the sum rule 



(5.7) 



4na 2 f3 J (log cosh ^) p®(s) dS = -1 (5.8) 

where p^ 2 \s) =< p(0)p(r) > is the charge pair correlation function. For a flat two- 
dimensional system (a — > oo), (|5.8|) reduces to the usual second Stillinger-Lovett 
sum rule 



r '(r) 27rr dr — — 1 



(5.9) 



6 DEBYE-HUCKEL APPROXIMATION 

Although the Debye-Hiickel approximation can be formulated for any many- 
component plasma, for the sake of simpler notation we consider only the case of a 
one-component plasma with number density n and particle charge q (the background 
has a charge density —qn). The Debye-Hiickel approximation is expected to be valid 
in the weak coupling limit f3q 2 <C 1. 

6.1 Correlations 

The average electric potential <p(s) at distance s from a given particle obeys the 
Poisson equation 



A0(s) = -2vrg \nh(s) + 5 {2 \s 



(6.1) 



while the pair correlation function h(s) approximately obeys the linearized Boltz- 
mann equation 
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h(s) = -(3q 4>(s) 

Thus, we obtain the usual Debye-Hiickel equation 



(6.2) 



(A-/t 2 )0(s) = -27iq5 (2 \s) (6.3) 
where the inverse Debye length k is defined by 

k 2 = 2ii(3nq 2 (6.4) 
and here the Laplacian is given by (|3.1|) . In terms of r = s/a, (|6.3|) becomes 

( —J— -^-sinhr-f- - K 2 a 2 ) = -2^ 2 \t) (6.5) 
ysmhr dr dr J 

The solution which vanishes at infinity is 



= qQu{coshr) 



(6.6) 



where Q u is a Legendre function of the second kind [ 1 1 and 



1 /l \ 1/2 
v= -- + (~ + ,V) (6.7) 

such that v{y + 1) = n 2 a 2 . The pair correlation function is given by ( |6.2| ) and the 
charge pair correlation function is 



q 2 n 2 h(s) = -i (^-j Q u (coshr) 



(6.8) 



At large distances s, (|6.8|) has the exponential decay — exp[— [y + l)s/a]. 



It can be checked that the screening sum rules (5.1) and Q5.8D are satisfied. 
Indeed, (|5.1|) results from |T0[ 



Q v [x)dx 



(6.9) 



As to (|5.8|) , one can use the Legendre equation (|6~5| ) for showing that [//(zz+^J^^a; 2 — 
l)dQ u (x)/dx is a primitive of Qy(x), and by successive integrations by parts of the 
left-hand side of ( |5.8| ) one reduces it to (|6.9j). 

Thus, like in the case of a flat system, the linearized Debye-Hiickel approximation 
exactly obeys the sum rules. This means that the system is a conductor. 
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6.2 Thermodynamics 

From the correlation function (|6.8|) , one can obtain all the thermodynamic 
functions, which are reasonably defined without boundary effects as follows. 

The excess internal energy per particle u exc (i.e. the potential energy per particle) 
is half the energy of a particle in the potential created by the other ones : 



Uexc — ~Q lim 
2 s^o 



[s) + glogtanh — 



(6.10) 

Using the behavior |K| of Q v (coshr) when z = coshr = cosh(s/a) is close to 1, 



Qu{z) ~ -\ log Z -^- - 7 - #/ + 1) (6.11) 



one finds 



Uexc = ~j [7 + ^ + 1)] (6.12) 

^From ( |6.12| ), where v is a function of (3 and n through (|6.7|) and ( |6.4|) , one can 
obtain the other thermodynamic quantities of interest as defined by (|4.7|) . However, 
the integral leading to / cannot be expressed in terms of known functions. 

6.3 Limiting cases 

It can be checked that, in the flat system limit a — > oo, the above results lead 



to the known behaviors [TB 



q 2 

h(x) ~ -f3q 2 K (Kr) , u exc — (7 + log na) 

q 2 ( 1\ q 2 



fexc ~ -— (7 + ^ogna - -J , p exc — n (6.13) 

The other extreme case is when the density n becomes small for a given value of 
the curvature radius a. The excess thermodynamic functions are analytical in the 
density. The excess internal energy (|6.12|) has its density expansion starting as 

u exc = -^-/3q A a 2 n + --- (6.14) 
6 
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in agreement with the weak-coupling limit (3q 2 C 1 of ( f4.10|) . The other weak- 
coupling limits are 



7T 3 

-,4„2„ 



fe XC = -— PqWn + ■ ■ ■ (6.15) 



PPezc = -^(Pq 2 )W + ■ ■ • (6.16) 



7 ONE-COMPONENT PLASMA AT (3q 2 = 2 

This is an exactly solvable model in a variety of geometries. On a pseudosphere, 
the pair correlation function has been recently obtained by a mapping on a field 
theory jnj . Here, all correlation functions will be derived by adapting the methods 
which have been used for the previously studied geometries |15[ . 



7.1 Correlation functions 

As explained in Section 4.3, one expects to find the same correlation functions 
either with the interaction — logtanh(s/2a) or with the interaction — logsinh(s/2a). 
It turns out that it is simpler to use the second one, as follows. 

For a one-component plasma with a — log sinh(s/2a) interaction, the potential 
energy can be obtained by adapting Q4.14p . Since 



\z — z' 1 



Alog^ — l - = 27i6 {2 \z,z') (7.1) 



the particle-background interaction energy arising from the first term of ( 4.14 ) is 
W(r) such that 

AW(r) = 2vrg 2 n (7.2) 
The solution of ( |7.2j ) is, up to an additive constant, 

W(r) = -2ixq 2 na 2 log |l - £^ (7.3) 
(choosing this W depending only on the distance r to the center of the Poincare disk 



is just a matter of convenience). The second term of (|4.14j) is another contribution 
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to the one-body potential. The particle-particle interaction comes from the first 
term of Q4.14Q . The background-background interaction is a constant. All together, 
the total potential energy is 



H = -q 2 ^og 

i>j 



I 



-£!(w + i)x>g(i-^) 



(7.4) 



up to some additive constant (actually an infinite one in the thermodynamic limit) 
which is irrelevant for the calculation of the correlations. 

The formalism for dealing with a potential energy of the form ( [7.4|) , in the canon- 
ical ensemble, has been previously developed [15]. When (3q 2 = 2, the Boltzmann 
factor is 



e-f> B = C 



2-irna 



ni'-'iii'-s, 



i>j 



where C is a constant. The Vandermonde determinant identity 



II (* 

i>j 



det 



1 


i 


1 


Z\ 


Z2 


Z3 








4 


z 2 


Z 2 
z 3 



allows to rewrite ( |7.5|) as the squared modulus of a Slater determinant 

det {^(^)}i, i= i |2 ,3,.. 



-PH 



c 



where 



/ 2 \27rna 2 + i 



(7.5) 



(7.6) 



(7.7) 



(71 



Thus, ( |7.7| ) has the form of the squared modulus of the wave function of a system of 
independent fermions occupying the mutually orthogonal one-particle wave functions 
ipj(z), and computing the n-body densities is a standard problem. The projector on 
the functional space spanned by the functions ipj is, in the thermodynamic limit, 



< z x \P\z 2 >= J2 



i/; j (z 1 )ip j (z2) 



PLS\M*)\ 2 *s 

In terms of this projector, the n-particle truncated densities are 



(7.9) 
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n(z) =< z\P\z > 

n$ ) (z 1 ,z 2 ) = -|< zi\P\z2 >| 2 (7.10) 

TlP (Z U Z 2 , ■ • • , Zn) = (~) n+1 J2 < Z h\ P \ Z h > < Z in\ P \ Z h > 

(i,i2—in) 

where the summation runs over all cycles (i, i 2 ■ ■ ■ i n ) built with {1, 2, • • • , n}. 
Using (|2.5|) and ( |7.8|) , one computes the normalization integral 



and the series (|7.9|) can be summed into 



_2 \ 27rncr + 4 , 2 \ lirna 1 



<z 1 \P\z 2 >-n _ x 4W+i y 7A2 > 

( 1 - S12, 



Using ( [7.12[ ) in ( |7. 10| ) gives the n-body densities. 

The one-body density is found to be position- independent, with the value n such 
that the particle charge density qn is the opposite of the background charge density, 
as expected. 

After simple manipulations using ( |2.4j ), the two-body truncated density ( [7.1 0| ) is 
found to be, in agreement with ref. |l 



n 2 h(s 12 ) = n?\ Zl , z 2 ) = -—— j , lW+2 (7-13) 

[cosh(si 2 /2a)J 

where S\ 2 is the geodesic distance between the two particles. This formula can be 
recovered from the analogous result on a sphere JT| of radius R by replacing R by 
ia. For large values of the distance S\ 2 , ( [7.131 ) has an exponential decay. 

It can be easily checked that the screening sum rules (|5.1| ) and (|5.8|) are satisfied 
by (|7.13| ) (for calculating the integrals, it is convenient to use cosh(s/2a) as the 
variable). Thus, the system is a conductor. 

In the Appendix, it is checked that using the interaction — logtanh(s/2a) gives 
the same correlation functions. 
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7.2 Magnetic analogy 

As in the case of a flat system, there is a magnetic analogy. For a quantum 
particle of unit charge living on the surface and submitted to a magnetic field B 
normal to the surface fl6fl , the wave functions associated with the infinitely de- 



generate lowest Landau level are of the form ( |7.8j ), in an appropriate gauge, with 
2irna 2 + (1 / 2) = Ba 2 . For a system of independent fermions filling the lowest Landau 
level, the two-body truncated density is just (|7.13|) . 

7.3 Thermodynamics 

The excess internal energy per particle is 

u exc =^J h{s)q 2 v{s)dS (7.14) 
where h(s) is given by (|7.13|) , v(s) by ( |3.4|) , and dS by ( f4.3|) . Using the variable 



t = tanh 2 (s/2a), after an integration by parts one obtains 



- f_ f 1 

Uexc , I 

4 Jo 



;i - t) 4W - lj — (7.i5) 



dt 



This is related to an integral representation [|l(| of the ip function such that 



Uexc . 



7 + ^(4vma 2 + l) (7.16) 
In the a — >• oo limit, one correctly recovers the flat system result 



J ~ 4 (7 + lo S 4 ™ 2 ) ( 7 - 17 ) 
The other extreme case is the small-n expansion (indeed, u exc is an analytic function 
of n around n = 0) 



7T 3 



Uexc = — —q a n + 
o 



which agrees with the result ( 4.10 ) from the virial expansion with (3q 2 = 2. 

Since u exc is known only at one temperature {(3q 2 = 2), it is not possible to use 
( |4.7| ) for obtaining the free energy per particle and the pressure. Only the beginning 
of the virial expansion is available through (|4.6|) and (|4.8|) : 
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exc 



—7ia 2 q 2 n + ■ ■ ■ 



(7.18) 



(3p = n — 7iPq 2 a 2 n 2 + ■ ■ ■ 



(7.19) 



8 TWO-COMPONENT PLASMA AT (3q 2 = 2 



This is a system of two species of particles with charges q and —q. When j3q 2 = 2, 
this model is exactly solvable in a variety of geometries. In the present geometry of 
a pseudosphere, the correlations can be calculated, in the grand-canonical ensemble, 
by adapting the previously used methods [0, [T7|] . 

8.1 Correlations 

The correlations can be computed by noting that ( |4.13| ) or (|4.14j ) map the system 
onto a flat system in the Poincare disk. The last term of ( |4 . 1 3| ) or ( [4.14| ) gives to 
the Boltzmann factor exp(—f3H) a multiplicative contribution [1 — (r,/2a) 2 ] for each 
particle. Furthermore, in the computation of a partition function, the area element 
dSi can be expressed in terms of the area element d 2 Zi as dSi = [1 — (ri/2a) 2 ]~ 2 d 2 Zi. 
Thus, the original system with a constant fugacity ( maps onto a flat system with 
a position-dependent fugacity ([1 — (r/2a) 2 ] -1 . 

For that flat system, it results from previous work |1^] that the n-body densities 
can be expressed in terms of Green functions G ei£2 (z\, z%) {£% is the sign ± of the 
particle at Zj). For instance, the one-body density in the flat disk for particles 
of sign e at z is 47rCa[l — (r/2a) 2 ] -1 G ££ (z,z). The corresponding density on the 
pseudosphere has a supplementary multiplicative factor d 2 z/dS and is n e = 4ir(a 





(8.1) 



the one-body densities are 
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n £ = 4ir(aG ££ (z, z) 



(8.2a) 



Similarly, the two-body truncated densities are 



n SK^i^2) = -e 1 e 2 (4:Tc(a) 2 \G £l£2 (z u z 2 ) 



(8.2b) 



and higher-order correlations are given by sums over cycles similar to the ones in 



( |7.10 ), the elements of the cycles now being 4ir(aG 



SiSj K Z iy Z j 



The Green functions G ££ >(z : z') form a 2 x 2 matrix G. It was shown [[T7j that 



when the potential energy is of the form (|4.14 ), G obeys the equation 



4n(a 



1 - (r/2a) 2 
where p is the flat Dirac operator 



G(z,z r ) 



(8.3) 



p= 2 



d z 
ds 



U) 



c(2) 

and I the unit 2x2 matrix ; 8W at is the usual Dirac distribution in the plane. From 
T|) and ( p.3|) , one finds 



(ip+AuCa) G{z,z')=I5 {2) ( 



z, z 



(8.5) 



where 



2 \3/2 



-1/2 



■6) 



is the Dirac operator on the pseudosphere, and now 5^ is the Dirac distribution on 
the pseudosphere, obeying (|3~^ ) . Thus G is the Green function of I/) + 4:7r(a, on the 
pseudosphere, and the densities have the simple expressions (8.2) in terms of G. 

The operator (l/4ira)p appeared in the formalism as the inverse of a matrix 
kernel 



/ 



M{z,z') 





2a 



2a 
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associated to the potential energy ( |4.14|) 



pM(z, z') 



Ana ' K ' ' ~ f! " 1 
If, instead of ( 4. 14j ) , one uses the potential energy ( 4.13Q which is like the one of a 

system with ideal conductor walls, (|8.7|) is replaced by |L8|, [19 



/ 



M(z, z') 



V 



AO? ~ zz' 
2a 



2a 



•9) 



4a 2 - 



However, ( |3.8| ) remains valid (because 4a 2 — zz! has no zero inside the Poincare 
disk), and one finds for G the same equation (|8.3| ). It can be shown that in both 
cases ( |4.13|) and ( |4.14j ), G(z, z') should go to zero at infinity, i.e. when z goes to the 
boundary circle of the Poincare disk. Therefore, the correlations are the same ones 
for a two-body interaction =pg 2 logsinh(s/2a) or =pg 2 logtanh(s/2a). 
One element of the matrix equation ( |3.3|) is 



2d- z G ++ (z,z') + 



4:7i(a 



G-+(z,z') = 



5.10) 



1 - (r/2a) 2 

By combining ( |8.10| ) with the other equation relating G ++ and one finds an 

equation for G ++ alone : 



-4a 



r^_\ (47rCa) 2 
4a 2 ) z 1 - (r/2a) 2 



G ++ (z,z') = 4<KCa6ff at (z,z') 



(8.11) 



Similar equations hold for G and Since n £ is expected to be position- 

independent and n^^(zi, z-i) to depend only on the geodesic distance s between 
points Z\ and Z2, for obtaining these quantities it is enough to compute G(z,z') in 
the simplest case z 1 = 0. One looks for a solution of ( 8.11 ) depending only on r and 
one uses the variable u = cosh 2 (s/2a) = [1 — (r/2a) 2 ] -1 . Then, ( |3.11| ) becomes the 
hypergeometric equation 



d 



d 



u(l - u)— - - u— + Un(a 
du 2 du 



2\2 



G A 



u > 1 



(8.12) 



with the boundary condition that G ++ behaves like — (a\og(u — 1) as r — > (u — > 1). 
Furthermore, G ++ should vanish at infinity (u — > oo). The solution of ( |S.12j ) which 



satisfies these conditions is 10 
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~ , r(4vrCa 2 )r(47rCa 2 + 1) ^„a / 9 9 , , , 

= C« 1 w s V 2 !\, -m ( F (4vrCa 2 , 47<a 2 + 1; 87<a 2 + 1; u" 1 

1 (87r(,cr + 1) v 

(8.13) 

Using ( |S.10D gives JLO 



G-+(J,0) 



S r(8vrCa 2 + l) V 7 
xF (4vrCa 2 + 1, 4?rCa 2 + 1; 8vrCa 2 + 1; u -1 



(8.14) 



One also finds G (-2,0) = G ++ (z,0) and G+_(z,0) = — G- + (z,0). By using the 

expansions of the hypergeometric functions with respect to 1 — iz -1 , one can check 
that these Green functions have the proper flat-system limits as a — > oo for a fixed 
value of m = 4ira(. 

Using these results in QS.ID and (8.2) gives the one- and two-body densities. 
However, G ££ (z, 0) has a logarithmic divergence as z — > 0, and one should use 
G ££ (a,0) in fl3.2aj) , where a is a small cutoff, rather than G ££ (0,0). This amounts 
to replacing the point particles by small hard disks of diameter a. The final results 
are 



n± = C(4vrCa 2 



2a 



-27 - ip(4ir(a 2 ) - ^(47rCa 2 + 1) + 2 log — 



(the expansion [[T(| of (|8.13|) with respect to 1 — u 1 has been used), 

. 2 [r(47rCa 2 + l)] 4 1 



nil s 



nL 2 f( S ) 



[r(87rCa 2 + l)] 2 [cosh(s/2a)] 16 < a2 + 2 
1 \ 1 2 

F ( 4vrCa 2 , 4vrCa 2 + 1; 8vrCa 2 + 1; 



1 



(cosh I-) 2 < 



(2)?7 ^ 



(2)T/ n _ ^2 [r(4vrCa 2 + l)] 4 [tanh(s/2a)] 



[r(8vrCa 2 + l)] 2 [cosh(s/2a)] 16 < a2 + 2 

2 



F 47rCa 2 + 1, 4vrCa 2 + 1; 87rCa 2 + 1; 



(8.15a) 



(8.15b) 



(8.15c) 



(cosh^)^ 

For the calculation of higher-order n-body densities, the Green functions G(z, z') 
with arbitrary z' are needed. They would be obtained by a conformal transformation 
such that maps onto z', in the same way as in the case of a sphere H. 
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At small distances s, ( |8.15b|) and ( |8.15c| ) behave like in the flat system case 



IT?]] , i.e. as —(logs) 2 and l/s 2 , respectively. At large distances, they both have an 
exponential decay as =F exp[— (8n(a 2 + l)s/a]. 

The strong-curvature limit a — > is also a low-density limit (a 2 — > 0. In that 
limit (|8.15b| ) and ( |B.15c| ) should behave like ^( 2 /3q 2 v ~ =F4C 2 exp(— s/a) and they 
do. 

8.2 Thermodynamics 

Since the total density n = 2n + is given as a function of the fugacity ( by 
( |3.15a|) , the pressure can be obtained by integrating n = (d((3p)/d(. However, 



the integration cannot be performed in terms of known functions for arbitrary (. 
Explicit results can be obtained in limiting cases only. 



The flat system results [|T7[] can be recovered as a — ► oo for a fixed value of 
m = 4ira(. 

In the opposite limiting case ( — > 0, for a given value of a, one obtains 



n = 2n± = 2( + ( 16na 2 log — ) ( 2 + • • • (8.16) 



a 



which gives 



(3p = 2(+ (8vra 2 log— ) C 2 + • • • (8.17) 



a 



2a N 

f3p = n- (27ra 2 log— ) n 2 + ■ ■ ■ (8.1? 



and therefore 



(X 

This is in agreement with a direct calculation of the second virial coefficient ( |4.2| ) 
with a lower cutoff on s at s = a. 

9 CONCLUSION 

The present paper has already been summarized in the Abstract. Here are some 
final remarks. 
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Perhaps a possible way of computing a free energy without boundary effects 
would be to start with a finite domain of suitable shape with periodic boundary 
conditions and to take its thermodynamic limit. However, we have not been able to 
use this seemingly difficult approach. 



It is likely that the Kosterlitz-Thouless transition |20[ of a two-component plasma 
to a dielectric phase at low temperature does not exist on a pseudosphere. Since the 
potential — logtanh(s/2a) now goes to zero at infinity, the original argument about 
a balance between energy and entropy here says that breaking a pair of oppositely 
charged particles always leads to a lower free energy because the entropy wins. The 
two-component plasma is expected to remain a conductor even at low temperature. [] 

When (3q 2 = 2, i.e. when the temperature is twice the Kosterlitz-Thouless tran- 
sition temperature for a flat system, the two-component plasma is a conductor even 
in a flat space, and a fortiori should be a conductor on a surface of constant nega- 
tive curvature. It would be of some interest to check that conjecture by explicitly 
showing that the modified Stillinger-Lovett sum rule ( |5.8|) is satisfied by the corre- 



lation functions ( |3.15b|) and (|8.15c|) . This is left as an open problem, because of the 



technical difficulty of computing the relevant integral. 

Another, more fundamental, open problem is how to compute the pressure of the 
one-component plasma at the temperature at which the pair correlation function 
( [7.1 3|) is known. For a flat system, one would use the virial relation 

p = f3- 1 n-—q 2 I ' h(r)^- r dS (9.1) 
4 J dr 



Is there a generalization of this relation on a pseudosphere 



? 



3 An opposite case when there is no Kosterlitz-Thouless phase transition is a two-component 
plasma on the surface of a cylinder of finite radius and infinite length |2]J] . At large separations, 
the interaction potential behaves like the distance, the system behaves like a one-dimensional 
two-component plasma and is always in a dielectric phase p2j. 
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A APPENDIX : MORE ON THE CORRELA- 
TION FUNCTIONS OF THE ONE-COMPO- 
NENT PLASMA 



The correlation functions of Section 7.1 can also be obtained when using the inte- 
raction — logtanh(s/2a), by adapting the grand canonical formalism used in Section 
8.1 for the two-component plasma. 

Eq. ( |3.3|) could have been obtained from the integral equation [18, pij 



G{z,z') + / M(z,z") 



C 



7 G(z", z')d 2 z" 



1 



■M(z,z') 



(a.i; 



1 - (r"/2a) 2 v ' ' " Ana 

where M is the matrix (|8.7|) . Indeed, applying p to both sides of ( |A.1| ) gives ( |3.3| ) 
when ( |8.8| ) is taken into account. 

Eq. ( |A.1|) can be adapted to the case of a one-component plasma by keeping 
only the ++ components of the matrices G and M (from now on, G ++ will be 
called G). Furthermore, the particle-background interaction ( |7.3| ) brings into the 
position-dependent fugacity a supplementary multiplicative factor 



1 - 



r y 
2a) 



(A.2) 



where C comes from the additive constant in the particle-background interaction. 
Here C is important, having the property C — > +oo in the thermodynamic limit. 
With these modifications, (|A.1|) becomes 



G(z,z')+(e 



c 



1 - (zz"/Aa 2 ) 



A-Kna? — 1 



G{z", z')d 2 z" 



Aa 2 J 



Ana [1 - (zz'/Aa 2 )} 
(A.3) 

From (|A.3|) , it results that G(z,z') is an analytic function of z. For simplicity, we 
shall consider only the case z' = 0. If one looks for a solution of circular symmetry, 
G(z, 0) is necessarily a constant, and ( |A.3| ) becomes 



G 



1 + 



c~ 



n 



1 

Ana 



(A.4) 
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The n-body densities involve 4na(e c G which, in the thermodynamic limit C — > oo, 
becomes the background density n, independent of the fugacity (. Including in 
( |S.2aj ) and (8.2b) the multiplicative factor ( |A.2j ) for each fugacity gives n for the 
particle density, and 



z,0) 



-n 



l-(- 

2a 



Anna 2 



(A.5) 



for the two-body truncated density in agreement with ( [7. 13| ) . 

The above calculation illustrates that, in a grand-canonical approach to the 
one-component plasma with a fixed background, the bulk properties are fugacity- 
independent. 
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